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Abstract 


If the Hamiltonian for the motion of a charged particle in 
a magnetic field has one or more cyclic coordinates it may 
often he viewed as representing the motion of a particle sub- 
ject to a potential V * The use of V provides qualitative 
insight about the motion even in some cases where a solution 
of the motion cannot be obtained. Several examples using 
this concept are reviewed and discussed. 
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There exist two types^of classical particle motion which are 
relatively easy to handle in Hamiltonian form. One is the motion 
in a potential field, where the structure of the potential function 
not only gives the form of the Hamiltonian hut also provides quali- 
tative information about the resulting motion, and the other is the 
motion of a charged particle in a magnetic field. The purpose of 
this, work is to review examples in which motions of the second type 
can be reduced to the first type, with various related benefits. 

Perhaps the best known example of such a reduction is the work 
of Carl St&rmer performed around the turn of the century 

A.nd dealing with the motion of a charged particle in a magnetic 
dipole field. If p is the dipole moment generating this field, 
an appropriate vector potential is 

A = p r 2 sin© 'f (l) 

where (r, 6 , f ) are spherical coordinates with 0 measured from 

* 

the dipole axis and where 'f is a unit vector in the ‘-r -direction. 

For time- independent fields, a Hamiltonian for a particle of mass m 

( 2 ) 

( = y 9 where m Q is the rest mass) and charge q is then 

H = (l/2m) £p 2 +. r” 2 p^ +■ (r sinQ )“ 2 (p^ - qp sin^e/r) 2 ] ( 2 ) 

There exist two constants of motion - the total energy E and 
also the component p^ of the canonical momentum, since H does 
not contain explicit dependence on . These two constants, 
however, are not sufficient to provide an analytical solution for 
the motion and in fact, such a solution in general cannot be found. 
Sturmer 1 s great contribution was in noticing that the constancy of 
reduces the Hamiltonian to that of motion in a 2-dimensional 
potential 

H = (l/2m)(p 2 + r~ 2 p 2 )+ V(r, 8 ;p f ) 


( 3 ) 
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where 

V(r, 0 ) = (p^ - qp sin 2 0/r) 2 /2m r 2 sin 2 0 (4) 

The structure of V in the (r, 0 ) plane - conveniently repre- 
sented hy contours of constant V - will determine the nature of 
the motion, From its algebraic form, V is clearly non- negative ; 
furthermore, if p ^ and qp have opposite signs, V does not 
vanish at any finite distance from the origin. On the other hand, 
if the signs are identical, then V also vanishes along the line 

r = q |i sin 2 0 / (5 ) 

which represents a field line of the dipole field if is held 
constant , 

In addition, V in all cases tends to zero as r 00 , In 

the case where V vanishes on the field line ( 5 ) * that field line 
must be separated from infinity by a region in which V has finite 
values, and it therefore represents a localized ’’potential hole” • 

If the particles in question have sufficiently low energy, they will 
be trapped in this hole and can never reach infinity. 

Thus, by qualitative arguments based on the concept of the 

potential, St&rmer deduced the existence of trapped orbits in the 

dipole field. He actually went further and derived the ranges in 

energy and in arrival direction at various points on earth for 

y hi ch any arr iv ing part ic le s w ould have to lie t r appe d . St 8rme r ' s 

aim, by the way, was the study not of trapped orbits but of orbits 

arriving from infinity. Specifically, he wanted to trace the motion 

of particles responsible for the polar aurora, which he believed 

to originate on the sun, and his theoryjLed to a result which quail - 

tatively seemed to agree with observations, namely that if such 

particles had relatively low energies they could arrive only near 

the magnetic poles of the earth. Near the geomagnetic equator, he 

showed (assuming the earth's field to be a dijoole field) all low- 

v are m ovi ng m ^ 

energy particles reaching the surface trapped orbits. 
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Later investigations showed StfJrmer's explanation of the aurora 
to be incorrect,, but his theory fits very well the observed arrival 
of low-energy cosmic ray particles. Satellite observations since 
I958 have also shown that trapped orbits above the earth T s atmos- 
phere contain a large particle population, constituting the well- 
known radiation belt. 

Stomer's results are easily extended to axi symmetrical poloidal 

(3) 

fields, of which potential fields are a special subclass, as 

are magnetic mirror configurations. Let the magnetic field be given 
by Euler potentials^ as 

B = 7«C(r,e) * ( 6) 

Then H again reduces to the form (jj) , with V given by 

V = (p^ - q ^ ) 2 /2 m sin ^9 (7) 

This shows that with a proper sign of p^ , a "guiding field line" 
exists on which 


** ® P * / <1 

On that line V vanishes and around it particles may be trapped. 

For further analysis of particle motion in a dipole field the reader 

(o) 

is reierred to the very thorough work of Dragt ; both there and 

(l) 

in the , work of St8rmer contours of V in the (r, 9 ) plane 

are shown, normalized in such a way that one set of contours serves 
for all energies. 
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ONE - DIMENSIONAL POTENTIALS 


The preceding example concerns a case where qualitative insight is 
obtained about a particle f s motion in a magnetic field, even though a 
complete solution of the motion is not obtainable, by converting it 
to motion in a potential* in two dimensions. If the problem can be 
reduced to motion in a one dimensional potential, a solution is generally 
possible, although its explicit derivation may be complicated. The 
value of introducing a potential V in such cases is that it allows one 
to come to qualitative conclusions by quick inspection - and in particular, 
it allows the deduction of periodic character of the solution, which leads 
one at once to look for adiabatic invariants of the motion. 


As an example, consider the motion of a particle near a neutral 

(5 ) 

sheet, investigated by Sonnerup , in a magnetic field B which main- 
tains a constant direction hut with a magnitude varying linearly with 
distance along an axis perpendicular to B » In cartesian coordinates 


B = B Q (x/l) y (8) 

where L is a scaling distance and where B vanishes and reverses its 
direction at the plane x = 0 .An appropriate vector potential is then 


A = - (x^ Bq/ 2L) \ 


(9) 


from which 

H = (l/2m) [p 2 + p 2 + (P z + <Tx 2 ) 2 ] (10) 


where 

= 1 B q / 2L 


cr 
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Clearly p and p are constants of the mot ion * The y- component of 
the motion may he completely separated and will therefore he omitted here, 
leaving 

H - p^/2m + P 7^ 


where 


V 



2\i/ 

(Tx )/2m 


( 12 ) 


This may he completely solved in terms of elliptic integrals (as Sonnerup 
has done), hut quick insight may he gained hy merely examining the form of 
y „ As before. V is non-negative, hut in this case it tends to infinity 
as x -* t , showing that V represents a potential well and that the 

motion is therefore always periodic a 

If p and <T have opposing signs, V vanishes at two points, namely 
r z 

at - • 



These will he minimum points of the curve of V against x - the only 
ones, it turns out - and that curve will resemble the one shown in 
Figure (l-a) * If <Tp z > 0 , only one minimum point exists 

and it is located at x = 0 , as shown in Figure (l-b) . 

In the first case, a mode of motion exists in which particles remain 
confined to one side of the neutral plane (i.e. the plane x = 0 , on 
which B vanishes ) » Such particles oscillate hack and forth in one of 
the "side pockets" of Figure (l-a) in a mode of motion resembling the 
well-known guiding center motion (Figure 2-a) - In addition, berth forms 
of V admit solutions which cross the plane x = 0 . The difference 
will he that in the second case (Figure l-b) the z component of the 


velocity 
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v = (p + <T x^)/m 
z z 


(15) 


maintains a constant sign, while in the first case (V as in Fig. 1-a) 
reverses sign at the points at which V - 0 • Thus in the second 
case the particle's advance in the z-direction is continuous, as shown 
in Figure (2-h) , while in the first case the particle's motion in 
that direction includes occasional retrograde motion, as shown in Figure 
(2-c) . 


GENERALIZATIONS AND ADIABATIC INVARIAM'S 


A more general, form which includes (8) as a special case is 



t \ * 

B 0 (x) y 


(ik) 


Such a field can he handled in much the same fashion, with the vector 
potential ‘ 

A = A(x) z ( 15) 


being derived by indefinite integration from 


A(x) 


J* B q (x) dx 


(16) 


The best-known case of this type is that of a constant field 
in which case A = B Q x and (12) is replaced by 


B oi 


> 


V = (p z + q B Q x) 2 /2m (l?) 

Regardless of the sign of p y V always has the profile of a parabolic 

z 

well and the Hamiltonian reduces to that of a harmonic oscillator. This, of 
course, gives the well-known guiding center motion of a particle in a 
uniform field. 
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The guiding center motion and the motion across a neutral plane are both 
periodic,, Classical mechanics then shows^^^ that the so-called 
action integrals derived for such motions 


§ 


P dx 
x 


(18) 


(integration over one period) remain adiabatic ally invariant in the presence 
of slov? perturbations o For the constant field B q y , the result is the 
familiar magnetic moment while for the motion in the field (8) a more 
complicated invariant has been derived by Sonnerup* 

J 

If a curl- free electric field - Is added, the solution may 

still be reducible to a- motion in a potentials in particular, if 

depends only on x , the Hamiltonian still has the form (ll) and 

only the function V appearing there is modified. If 4> = ^ (x, y) 

the potential becomes two-dimensional, while if dependence on z is 

also allowed, p is no longer constant and all the advantage of 
z 

introducing V ' ' Is lost . 

An exception to this is a constant electric field E , for which 




= xE + y E f zE 
x J y z 


( 19 ) 


In that case, the motion in y still separates, and if the magnetic 

field is constant, the problem may still be reduced to motion in a 

one -dimensional potential. To do so one employs a canonical transformation 

to new variables (P., Q. ) , introduced by Gardner if B is given as 

a 11 “ 

Bq y , the generating function of the transformation is 


F(P, r) = B-P, x +- P y -t- P, z - FjP, / q 


( 20 ) 


A few additional steps reduce H to a form similar to that derived for a 
constant field, but the only net result is the addition of the guiding 

•f i 

center drift due to E „ 
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TWO - DIMENSIONAL FIELDS 


Consider the general class of fields independent, of z and lacking 
a z-component 

B - x B (x, y) + y B (x, y) (Zl) 

“ x — y 

(B and B must be related by zero- divergence condition). Such fields 

* y M 

can be expressed in terms of Euler potentials 

B = VCJC(x -> y)xV'z (22) 


where &£ is the indefinite integral 


(x, y) 



y) dx 




B i (y) 


and where is adjusted to satisfy 


^ oC /T>y = B x (x, y) 


(23) 


A field line of this field is then characterized by the constancy of 
(x, y) and of z, and the Hamiltonian is 

H = tl/2m)[p x 2 + p y 2 ] +• V(p 2 , X, y) (24) 

where p is a constant of the motion and 

n 2 

= j^P z ** 1 «<(x 3 y) J /2m 


V 


(25) 
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The motion thus reduces to a two-dimensional motion in a. non-negative 
potential in which all equipotentials follow magnetic field lines , with 
the third coordinate arbitrary and constant* 

magnetic 

As an example, conside the X-type^null line (or neutral line - a line 
along which B vanishes) obtained by adding an extra term ,to equation (8) 


B = B. (x/L) x i" B (y/L) £ (26) 

— JL d 

The appropriate Euler potential will satisfy 


q o( (x, y) 




(27) 


where C T = q B./2L » The field lines are a family of hyperbolas 

i 1 

sharing the same asymptotes (Figure 3) and the Hamiltonian is 


H = (l/2m) 




(28) 


All equipotential lines are thus field lines (and vice versa) but unless 

p is given, one cannot determine which of these lines represent high 
z 

values of V and which ones represent low values. For any given p^ , 
there will exist two conjugate hyperbolas on which V vanishes 


P 


z 



and these represent the bottoms of two valleys, separated by a pass at 
the origin* A particle starting from one of the valleys (broken line 
in Figure 3) needs a certain minimum of energy to climb out of it and 
to cross the pass. However, even when the energy exceeds this minimum. 



11 - 


the particle might fail to cross the pass: whether it. succeeds in doing 
so depends on the exact initial conditions a The motion may he investigated 
qualitatively and numerically^ hut a complete analytical solution is 
not available „ 

Other two-dimensional fields with some practical interest are those 
which can serve as models for the central portion of the earth s 
"magnetic tail" ; such fields differ from the one in equation (8) by 
small additional terms. For instance 


B 


b q (x/l) y + 


B l - 


( 50 ) 


has parabolical field lines and represents one such model, with B 1 

accounting for the average northward magnetic flux observed in the 

geomagnetic tail, due to the magnetic polarity of the earth's main 

field. This field has the interesting property that first order guiding 

( 12 ) 

center drifts in it (curvature and gradient) cancel identically 
when averaged over a single longitudonal excursion oi the particle. 

(] 3 ) 

Another model contains a string of alternating X-type and 

such behavior 

0-type null lines and^may be represented by 


B = B q (x/l) y + B^ x sin wy (51) 

Both these fields can be handled in the same way as the field of 
figure (5), but details are left to the reader. 


BOUNDARY FIELDS 

Grad ^ ^ ^ Investigated the transition between a field- free plasma and 
a plasma- free uniform magnetic field. The magnetic field in his example 
was 

B = B(x) z (5?) 
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with B(x) i 0 for any finite x and with 


B(o«o ) -> B q = const. 

B(- oe> ) — * 0 


The plasma density vas assumed to depend only on x , tending to a 
constant limit as x -*• - and to zero as «>«* . 

The vector potential for (32) can he written as 


A - 


A(x) 9 


(33) 


with A(x) derived by integration 


A(x) 


x 



(34) 


Because B j- 0, the sign of B (assumed to be positive) never changes , 
making a(x) monotonic ally increasing. The Hamiltonian of the motion in 
the (x, y) plane (the z-component of the motion may be separated) is then 


H 


(l/2m) £ 



¥ 




(35) 


and has a form similar to that of (ll), since is conserved. The 
potential V now has the form 


V = (p - 1 A ) 2 /2m (36) 

and the qualitative analysis resembles that of (ll). If p and qA 
have the same algebraic sign, then V increases monotonically with x , 
so that any particle for which this is true will penetrate to a maximal 
value of x and then roll back all the way to - 00 . 
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If the signs differ, there will exist for the given value of p an 
absolute minimum of V at which V vanishes, and particles of suffici- 
ently low energy will be trapped In the vicinity of that minimum- In 
conventional terms, orbits extending to x = - correspond to parti- 

cles arriving at the boundary of the magnetic field from the field-free 
region and being reflected at that boundary, while trappedjorbits repre- 
sent guiding center motion inside the region of appreciable magnetic 
field* 


The Hamiltonian formulation is especially useful in this case 
because - through the Liouville equation - it leads directly to the 
Vlasov equation for a plasma located in the given magnetic field, and 

( 15 ) 

this has been investigated by various workers * The method may be 
expended to fields having the form 

B = B (x) y + B (x) z (37 ) 

y “ z 

The vector potential can then be viewed as the sum of vector 
potentials derived separately for each component, the derivation in 
each case following the one given in equations (33) “ (3^-)« Writing 

A = A (x) y + A (x) z ( 38 ) 

- y ** z - 

one can derive the components by integrating the relations 

dA /dx = B dA /dx = - B (39) 

y 1 z z' y 

The Hamiltonian is still that of motion along the x axis with a poten- 
tial V , but now 


t 2 2 1 

(p - q A (x) ) + (p - q A (x) ) /2m (40) 

y y z z J ' 

and both p and are constants of the motion* Again, the^transi- 

y 2 

tion to treatment of the Vlasov equation is straightforward: the general 
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solufcion of that equation will he a function of the constants of motion 

(-n -p e) and must he constructed in a way which makes the density of 
'•ky* 

electric charge and current agree with the postulated fields and avoids 
ne g at ive part i c le dens it ie s <. 


The important case where B represents a monotonic ^transition 
between two constant fields has "been analysed hy Alpers ° Such a 

transition may serve as a good model, for the magnetopause - the ^harp 
transition observed between the geomagnetic field and the interplanetary 
plasma - and also for "tangential discontinuities" observed by spacecraft 
in interplanetary space ^ 0 In the former case an electric field - V$(x) 
probably exists at least on one side of the boundary, ciue to tangential 
flow of the plasma there , and this adds one extra term to (40 ) • 


For large values of x , in these examples , the field tends to 
constant values and by (59 ) > i- n this limit 


A (x) 

y 



x B 


zl 


A (x) 
z 


A 


zl 


x B 


y! 


( 41 ) 


A similar limit exists as x-*- - °° , with index 2 replacing index 1 * 

In both cases, because of the factor x , the contributions of A to V 
■will be relatively large * Of course, the values of p and p^ could he 
of comparable magnitude, causing V to have its minimum at some large value 
of x and leading to guiding center motion in the vicinity of that x • 

For most values of p and p , however, when x is large the terms 

y z 

proportional to it in (^d r ^lLominate V and cause potential walls 

to rise parabolically as x ± , leaving between them a "valley" 

with a floor which may be quite irregular. Particles of sufficient 
energy will bounce back and forth between such walls, thus crossing and 
recrossing the boundary region in a manner somewhat similar to that of 
particles near Sonnerup’s neutral sheet (note however that nere the motion 
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is not confined to a single plane). Such particles trapped in 
the vicinity of the magnetopause seem to he responsible for the current 
density required to maintain it. 

Another example of a field having a form for which the motion of a 
particle is readily reduced to motion in a potential V Is the rotating 
field 


B = B ( £ sin -f- z cos >x ) (4 2) 


The details are left to the reader. The motion of charged particles 
can also he solved in a somewhat more general field which in addition 
to (42) also contains a constant magnetic field in the x direction 
and that case, too, reduces to motion in a one -dimensional potential, 
hut the potential is in velocity space and its derivation is not as 
straightforward as in the examples given here. Details may he found in 
the work of Lutomirski and Sudan 

CONCLUSIONS 

The preceding examples show that if the Hamiltonian for a charged 
particle in a magnetic field contains one or two cyclic variables, it can 
often be cast in a form representing motion in a one -dimensional or a 
two-dimensional potential V • The structure of V then often provides 
qualitative insight into the range of motions which can exist. If the 
motion reduces to a one -dimensional case, the existence of potential 
wells suggests periodicity, allowing adiabatic invariants to be deduced. 
In two-dimensional motions the derivation of complete solutions may 
prove impossible but qualitative conclusions might still be deduced from 
V , as in St firmer 1 s calculation. 

An interesting - and unanswered - question is how valid are the 
qualitative properties deduced from V when the strict conditions which 
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allow it to be introduced (usually some symmetry) are relaxed. In the 
one -dimensional case, adiabatic invariants deduced from V are expected 
to be approximately preserved even in the presence of perturbations. The 
derivation of adiabatic invariants, however, rests on classical pertur- 
bation theory and requires the availability of the solution for the 
unperturbed motion. In two-dimensional cases - like St&rmer's theory - 
such a solution is not available, and it is not clear how the qualitative 
theory can be extended to perturbed cases. Observations suggest that 
cosmic ray particles in the geomagnetic field conform quite well to 
Stfrrmer's theory, in spite of an appreciable asymmetry of the field: this 
suggests that an extension of that theory might well exist, although 
none is known at the present time. 
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CAPTIONS TO FIGURES 


Figure 1 


Potential veils for the motion of a particle in a magnetic 
field B q (x/l) y „ 


Figure 2 


Possible modes of motion for a particle in a mag netic field 
b q (x/l) £ o 


Figure j5 


Potential contours for the motion of a particle near a 
magnetic null line,, 




~£ 0 - 



Figure 2- a 



Figure 2-i) 



Figure 2-c 
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Figure 3 


